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Abstract- By using the data for the known geopotential models by means 
of artificial satellite, the central moments of inertia of the Earth are determined. 
For this purpose, it was used the value H = 0.00327369±9.8-10~* for dynamical 
flattening of the Earth [7] . The results obtained indicate that the pole of inertia 
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1 Introduction 

The first artificial satellite of the Earth, Sputnik- 1, was launched on 4 October 
1957. Studying the trajectory of the following satellites (Sputnik-2 and Sputnik- 
3), D. King-Hele has determined the zonal coefficient J2. The value J2 = 1.084- 
10 ~^ was quite close to that calculated by terrestrial measurements. The US 
satellite, Vanguard-1, launched in March 1958, made it possible for the first 
time the assessment of discrepancy between the ellipsoid and geoid. The value 
of J4 was obtained in the same year and the first odd zonal term in 1959 by Y. 
Kozai. In 1961, W. Kaula produced a complete model of degree 4, involving all 
the coefficients Cim and Sim of the associated Lagrange function P22 • From this 
moment, the information about the gravitational field of the Earth has become 
more numerous and accurate. 

The first data from satellites have been used in the development of geopoten- 
tial models from the early 1970. The SAO-SE model (Smithsonian Astrophys- 
ical Observatory - Standard Earth), established in 1966, used in 1972 the first 
laser-ranging measurements to establish satellite distances. The GEM model 
(Goddard Earth Model) was established by NASA's GSFC (Goddard Space 
Flight Center) in the United States as a reaction to the classified US military 
models. The first model, GEM-1, was published in 1972, expanding the poten- 
tial to degree 12. He then followed the whole series of geopotential models to the 
GEM-10 (developed up to order 20). Subsequently it was developed the model 
EGM, as a resuh of the coUaboration between GSFC-NASA, NIMA (National 
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Imagery and Mapping Agency) and OSU (Ohio State University). In 1996 came 
EGM96S, of degree and order 70, with data provided solely by satellites, and 
EGM96, of degree and order 360, adjoining geophysical data |T|. 

In this paper, using the data provided by the SE-2 geopotential models from 
SE series, GEM-5 to GEM-10 models from GEM series and the EGM96 model, 
the Earth's moments of inertia are calculated. Although at most six digits 
are accurate, to compare the results more easily with the values obtained by 
Erzhanov and Kalybaev [3], the calculations are performed with nine digits. In 
order to determine the polar moment C" it was used the equation obtained by 
Prof. leronim Mihaila W- In fact, for all the nine models used, the value of C" 
in Tables 4 and 6 coincides with the value of C" obtained by considering 



is obtained from the constant of precession. We used here for dynamical flat- 
tening the value H = 0.00327369 ± 9.8 • lO'^ [7]. Thus, it demonstrates that 
the choice made in [5, namely H' = H, is valid until the order of 10~^. The 
orientation of the ellipsoid of inertia is described in Table 7. The data from 
satellites on the gravitational potential indicate that the equatorial principal 
moments of inertia of the Earth are not equal (see Table 6) and it also shows 
us that the polar axis does not coincide with the axis of rotation. The pole of 
inertia Pi remains near the Conventional International Origin. 

2 Representations of the geopotential 

In the theory of the movement of the Earth's artificial satellites it is chosen as 
a reference system the geocentric system O^jyC, the axis OC of the system being 
given by the position of the Conventional International Origin. The origin plan 
for longitude, 0£,(, is the plan of the Greenwich meridian. 

In polar coordinates, the expression of the geopotential is [3] 
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where 

n 

Yn{(p, A) = ^ Pli"^\cosip)[AnmCosmX + BnmsinmX]. (1) 

m— 

Here, A is longitude and Lp is the geocentric latitude. The symbol V indicates 
that the integration should be extended to the whole volume of the Earth. The 
coefficients and -B„m from ([T]) are 
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while P„ and P,!™'' are respectively the conventional zonal harmonics of n*'' 
degree and the associated function of Legendre of n*'* degree and m*'' order. 
We mention that the recommended geopotential form by U.A.I. [5] is 
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where the harmonics coefficients of the geopotential are 
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There are three kinds of Legendrians in use: the conventional Legendrian 

"™ ~ (n+m)! ' 



when hnm — ^ the normalized Legendrian when hnm = 1 and the fully 
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normalized Legendrian when hnm = (2n + [H]- We use the conventional 
Legendrian. 

The connection between the harmonics coefficients of the geopotential from 
([3]) and the coefficients ^ is given by the relations 

~ Mae" 

where M is mass of the Earth. 

Often, the geopotential is defined by the next expression: 
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where the relations between the coefficients /,„, Inm and the constants A„m with 
the coefficients ^ are (see [3]) 
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then the series ([3]) becomes (see 
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If we note the following terms 
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then the polynomials Pm and Pnm may be written as follows 
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The geopotential models are characterized by some constant values (see Ta- 
ble 1), called the universal constants of geopotential. They include the equato- 
rial radius of the Earth (ce), the geocentric gravitational constant (GM) and 
the geometrical flattening of the Earth (/e). 

In the following, it is necessary to know the dynamical flattening of the Earth 
(H) . We use the value of H calculated from the constant of precession |7j 

H = 0.00327369 ± 9.8 • 10"^ (5) 



MODEL 


ae[m] 


fe 


GM ■ 10-^'^[m^ ■ s-'^] 


SE-2 
GEM-5, GEM-6 
GEM-7, GEM-8 
GEM-9, GEM-10 
EGM96 


6378155.0 
6378155 
6378137.8 
6378139.1 
6378136.3 


1/298.255 
1/298.255 
1/298.7925 
1/299.7925 
1/298.257 


3.986013 
3.986013 
3.9860013 
3.9860064 ±0.02 
3.986004415 



Table 1. Universal constants of geopotential 
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3 Determination of the moments of inertia of 
the Earth 



Using the expressions ^ of the harmonics coefficients of the geopotential, one 
finds the relations between the first coefficients and the moments of inertia of 
the Earth A', B', C , D' , E' , F' in the system Oir^Q, namely 

(6) 
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For the moments of inertia of the Earth, the following notations were used: 

A' = I p{i,flX){rf +e)dv, 
Jv 

B' = [ pii,V,Oie+C')dv, (7) 
Jv 

Jv 

D' = I pitvXKvdv, E' ^ I p{^,ijXK^dv, F' ^ I p{^,ijX)^vdv, 



where p is the density. On the other hand, from (|4|) and ([6|) the following 
relations between the coefficients In, Cnm, Snm and the moments of inertia of 
the Earth are obtained: 

4MalC22 = B'-A', 

MalC2i - E', (8) 

MalS2i = D', 
2MalS22 = F'. 

The five relations ([S]) are insufficient to determine the six moments of inertia ([7]) 
of the Earth. The system would be complete if another independent equation is 
added. 

Erzhanov and Kalybaev [3] had the idea to use the following expression for 
the sixth equation of the system 

H' = ^[2C'-{A' + B')], 
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taking the H' = H without a motivation for this approximation. To avoid it, 
Prof. I. Mihaila deducted an equation for calculating the polar moment C", 
using for this purpose the expression of -ff [7]. 
With this equation of the polar moment: 

ax^ + bx'^ + cx + d ^ 0, (9) 

where the coefficients a, b, c and d are respectively 

a = 8i^^ 

b = 8H^a', 

c = -2H{3-4H)a'^ + 2H{3-2Hf{ D'^ - E'^ - F'^), 

d = -2(1 - H)a'^ + (3 - 2Hfa' C'^ ~ ^'^ - D'^ ~ E'^ - F'^) + 

+ (3 - 2Hf{2D'E'F' + ^^S'^ + ^-—^D'^), 

a' = -2J2Mal, 

b' = B' - A' ^ 4MalC22, 

the system of six independent algebraic equations ([5]) and © for the six searched 
moments of inertia A', B' , C", D' , E', F' is obtained. 

If the values of harmonic coefficients of order n = 2 for the geopotential and 
the dynamical flattening of the Earth are known, then the normalized moments 
of inertia A , B , C , D , E , F can be determined easily. We note here A = 
A' /Mai, etc. The SE -2 model is completed by the SE - 2' model, where the 
harmonics coefficients Cii = -0.001196 • 10^^ ^nd ^21 = -0.003466 • 10"^ ^^^^ 
calculated by Erzhanov and Kalybaev j^. 



MODEL 


C20 ■ 106 


C21 • 106 


521 • 106 


C22 • 10^ 


522 • 10^ 


SE~2 


-484.16596 


0.00000 


0.00000 


2.41290 


-1.36410 


SE~2' 


-484.16596 


-0.001196 


-0.003466 


2.41290 


-1.36410 


GEM - 5 


-484.16620 


-0.00120 


-0.00870 


2.42820 


-1.36020 


GEM - 6 


-484.16610 


-0.00090 


-0.00120 


2.42510 


-1.38830 


GEM -7 


-484.16460 


-0.00310 


-0.00090 


2.43030 


-1.39460 


GEM -8 


-484.16460 


-0.00010 


-0.00030 


2.43450 


-1.39530 


GEM - 9 


-484.16555 


-0.00021 


-0.00406 


2.43400 


-1.39786 


GEM - 10 


-484.16544 


-0.00104 


-0.00243 


2.43404 


-1.39907 


EGMm 


-484.16537 


-0.000187 


0.001195 


2.43914 


-1.40017 



Table 2. Normalized harmonics coef ficients of the geopotential (see[3], [4]) 
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MODEL 


^' 


b' 


c' 


SE-2 


0.329619974 


0.329626204 


0.330705717 


SE-2' 


0.329619974 


0.329626204 


0.330705717 


GEM-5 


0.329620259 


0.329626529 


0.330706023 


GEM-6 


0.329620507 


0.329625671 


0.330705717 


GEM-7 


0.329619038 


0.329625314 


0.330704801 


GEM-8 


0.329619033 


0.329625319 


0.330704801 


GEM-9 


0.329619643 


0.329625927 


0.330705412 


GEM-10 


0.329619643 


0.329625927 


0.330705412 


EGM96 


0.329619636 


0.329625934 


0.330705412 



Tables 3. a. The normalized moments of inertia of the Earth 



MODEL 


13' ■ 10*5 


'e' ■ 10*5 


f' ■ 10^ 


SE-2 








-1.761045528 


SE-2' 


-0.004474587 


-0.001544029 


-1.761045528 


GEM-5 


-0.011231652 


-0.001549193 


-1.756010649 


GEM-6 


-0.001549193 


-0.001161895 


-1.792287593 


GEM-7 


-0.001161895 


-0.004002083 


-1.800420858 


GEM-8 


-0.000387298 


-0.000129099 


-1.801324554 


GEM-9 


-0.005241437 


-0.000271109 


-1.804629500 


GEM-10 


-0.003137117 


-0.001342634 


-1.806191603 


EGM96 


0.001543100 


-0.000241400 


-1.807607613 



Tables 3.b. The normalized moments of inertia of the Earth 

In Tables 3, using the models of geopotential from Table 2, we evaluate these 
normalized moments. 

Further, if the mass M and the equatorial radius of the Earth arc known, 
then one can determine the central moments of inertia A' , B', C , D' , E' , F' . 
In Table 4. a. and Table 4.b., using the data from Table 1 and the value of 
the gravitational constant given by the lAU (1976) System of Astronomical 
Constants, namely G = 6.672 • 10~^^m'^/c(7~^s~^, these moments are evaluated. 

Once the values of the moments ^ known, the principal moments of inertia 
A, B, C can be determined by solving the secular equation (see [5], [5]) 



A(9) 



A' -q 
-F' 
^E' 



-F' -E' 
B' -q -D' 
-D' C -q 



= 



(10) 



The roots gi, q2, q^ of equation 
A, B, respectively C. 



10|) represent the principal moments of inertia 
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MODEL 


A' 


■ lQ-''^'^[kg ■ m^] 


B' ■ lQ-'^^'^[kg-m^] 


C" • 10-^qfc5 • TO^] 


SE-2 




8.010992630 


8.011144042 


8.037380227 


SE-2' 




8.010992630 


8.011144042 


8.037380227 


GEM-5 




8.010999557 


8.011151941 


8.037387664 


GEM-6 




8.011005584 


8.011131088 


8.037380227 


GEM-7 




8.010903161 


8.011055690 


8.037291025 


GEM-8 




8.010903040 


8.011055812 


8.037291025 


GEM-9 




8.010931380 


8.011084104 


8.037319434 


GEM-10 




8.010931380 


8.011084104 


8.037319434 


EGM96 




8.010920187 


8.011073251 


8.037308375 


Tables 4. a. The moments of inertia of the Earth 


MODEL 


D' 


• 10^^"[fcg-TO^] 


E' ■ lQ-''^'^[kg ■ m^] 


F' ■ 10-^^[/cg-TO^] 


SE-2 










-4.279996317 


SE-2' 




-1.087491183 


-3.752566228 


-4.279996317 


GEM-5 




-2.729709620 


-3.765116622 


-4.267759686 


GEM-6 




-0.376511662 


-2.823837454 


-4.355926167 


GEM-7 




-0.282381394 


-9.726470762 


-4.375656585 


GEM-8 




-0.094127316 


-0.313757105 


-4.377852885 


GEM-9 




-1.273855995 


-0.658891032 


-4.385892467 


GEM-10 




-0.762431051 


-3.263079405 


-4.389688933 


EGM96 




-0.375027747 


-0.586687176 


-4.393124297 



Table 4.b. The moments of inertia of the Earth 



MODEL 


A 


B 


C 


SE-2 


0.329619513 


0.329626665 


0.330705717 


SE-2' 


0.329619513 


0.329626665 


0.330705717 


GEM-5 


0.329619801 


0.329626987 


0.330706023 


GEM-6 


0.329619939 


0.329626239 


0.330705717 


GEM-7 


0.329618555 


0.329625797 


0.330704801 


GEM-8 


0.329618549 


0.329625803 


0.330704801 


GEM-9 


0.329619167 


0.329626403 


0.330705412 


GEM-10 


0.329619160 


0.329626410 


0.330705412 


EGM96 


0.329619148 


0.329626422 


0.330705412 



Table 5. The normalized principal moments of inertia of the Earth 

The values of the Earth's normalized principal moments of inertia A, £?, 
C, where A = etc., obtained by the solving of the equation (jlOp with the 
geopotential models from Tables 3 are presented in Table 5. In Table 6, using 
the data from Table 5 and Table 1, are evaluated the principal moments of 
inertia A, B, C. 
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MODEL 


A ■ lQ-''"''^[kg ■ m'^] 


B ■ IQ-'-^^kg ■ m^] 


C ■ 10-^^[fc.g-m^] 


SE-2 


8.010981426 


8.011155246 


8.037380227 


SE-2' 


8.010981426 


8.011155246 


8.037380227 


GEM-5 


8.010988425 


8.011163072 


8.037387664 


GEM-6 


8.010991779 


8.011144893 


8.037380227 


GEM-7 


8.010891423 


8.011067429 


8.037291025 


GEM-8 


8.010891277 


8.011067575 


8.037291025 


GEM-9 


8.010919812 


8.011095672 


8.037319434 


GEM-10 


8.010919642 


8.011095842 


8.037319434 


EGM96 


8.010908325 


8.011085109 


8.037308375 



Table 6. The principal moments of inertia of the Earth 

As seen from Tables 3 and 5 or from Tables 4 and 6, it is noticed that C" 
coincides with C. For all the nine geopotential models used, the value of C 
found here coincides with the value of C" obtained by considering 

H'^^['2C'~{A' + B')] = H 

as in the work [3^. It is thus demonstrated that the choice made by Erzhanov 
and Kalybaev, namely H' = H, is valid until 10~^. 



4 Orientation of the ellipsoid of inertia 

Let Oxyz be the system of the Earth's principal axes of inertia, whose coordi- 
nate axes are chosen so that 



Jv 

B ~ p{x,y, z){z^ + x'^)dv, 
Jv 

C = / P{x,y,z){x'^ +y^)dv, 
Jv 

p{x,y,z)xydv — / p{x,y, z)xzdv — / p{x,y,z)yzdv^O. 



The orientation of the system with respect to OS,riC may be given by the Euler 
angles. We use the notations from [3j 

P={O^N),a^{ONrOx)n={OC?dz), (11) 

where ON is the intersection between the plans OS^rj and Oxy, called the line 
of nodes. 
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Let {p(,P7j,Pc), {q$,qn,<lc)' (''4 1 i ''C ) direction cosines of the axes 

Ox, Oy, Oz in respect with O^-qC,. They are the projections of the unit vectors 
p, q, r of the principal axes in the system 0^r](. On the other hand, the 0^r]( 
system overlaps Oxyz by three rotations Rp, Rj, R^. The direction cosine have 
the foUowing expressions: 





— cos{x,£,) 


= cos/Scosa — 


sinPsinacosj, 




= cos{x,ri) 


= sinf^cosa + 


cosPsinacos^, 


PC 


= cos{x,C,) 


— sinasiw-f, 




9? 


= cos(y,^) 


= —cosf3sina 


— sinPcosacos^, 




= cos{y,i]) 


— —sin/3sina 


+ cosPcosacos^, 


ic 


= cos{y,C) 


~ cosasinj, 






= cos{z,^) 


— sinPsinj, 






— cos{z,r]) 


— —cosf3sin"f, 






= cos(z,C) 


— COSJ 





(12) 



and satisfy the orthogonahty conditions 



pp = qq = rr = l, (13) 
pq = pr = qr = 0. 



The direction cosines are given by the relations 

111 li2 la 1 



5ii 5a Sa ^6f^ +Sh + Sl 



(14) 



where Sii,5i2, 6i3, with i=l,2,3, are the cofactors of the elements in row i of the 
determinant A which appears in the secular equation (|10p. q being successively 
replaced by the principal moments of inertia A, B and respectively C (see [5], 
[S]). In the relation (HH), we have (711,712,713) = {Pi,Pv^Pc)' (721,722,723) = 
(9?,'?r,,<7c)' (731,732,733) = {r^,rr„r^). 

In determining the orientation of the principal axes of inertia, besides the 
orthogonality conditions ([13]), it is also required to fulfill the condition that the 
system Oxyz to have the same orientation as the system 0£,r](, namely (see [2]) 



Pi Pn PC 

qr, qc 

rc_ 



(15) 



Once determined the direction cosines, it can be obtained the three Euler's 
angles by the relations ([T^ . Thus, for determining the angle 7 it is used the 
following relation 

COS7 — rQ, 
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and to obtain the angle a, the following relations are used 

Pq ~ sinasin"/, 
qQ = cosasinj. 

For the angle /3 we have the following formulae 

sinfisin-i = p^g^ - g^p^, 
cos/Ssinj = p^q^ — q^P(^- 

If 7 = 0, then the axes Oz and 0^ coincide and the plan Oxy coincides with 
0^77. The problem of the orientation for the system Oxyz is reduced in this 
case to the problem of the orientation for the plan system Oxy in relation to 
Oir]. We have from p^ : 

= cosj3cosa — sinPsina = cos{l3 + a) = cosX, 

Pfj — sinj3cosa + cosjSsina = sin(/3 + a) ~ sin\, (16) 

= —cosl3sina — sinf3cosa ~ ~sin{j3 + a) = —sinX^ 

q^i = —sin/3sina + cosjScosa — cos{j3 + a) — cosX, 

where it was noted by A the angle between axes Ox and O^. The expressions 
give us the ellipse of inertia orientation in the plan 0^77. 

The orientation of the ellipsoid of inertia corresponding to geopotential mod- 
els is given in Table 7. a and Table 7.b. From the model SE — 2 with 7 = was 
determined the orientation of Oxy in relation to O^rj. This orientation is used 
as a standard for the other models. 

In Table 7.b, the longitudes and the geocentric latitudes of the ellipsoid of 
inertia axes are also given. For the determination of the coordinates A^i and ipA 
of Ox, we have the following relations 

p^ = cosXacoslpa, 

Pri = sinXACOSlfiA, 

Pq — sirup A- 

Similarly, for the determination of the coordinates Xb and (ys^ of Oy, we have 
the following relations 

q^ = cosXbcos(Pb, 
qri = sinXBCOSifB, 
= sirup B- 

For Oz, the longitude Ac is determined from the following formula 
Ac = 360° - (90° +^), 

where /3 is given by (|lip . 
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MODEL 


a° 


P° 


7 (in arcsec (")) 


SE-2 









SE-2' 


4.3 


-19 


0.9 


GEM-5 


-7 


-7.9 


2.2 


GEM-6 


19.4 


-36.8 


0.4 


GEM-7 


58.7 


-73.6 


0.8 


GEM-8 


3.5 


-18.4 


0.1 


GEM-9 


-11.9 


-3.0 


1.0 


GEM-10 


8.2 


-23.1 


0.7 


EGM96 


23.7 


-38.6 


0.3 



Table 7. a. The orientation of the ellipsoid of inertia 



MODEL 


Xa° 




Ab° 




Xc° 


(2 - fc) 


SE-2 


-14.7 





75.3 













SE-2' 


-14.7 


1.9- 10-5 


75.3 


2.5-10- 


-4 


251 


0.9 


GEM-5 


-14.9 


-7.3-10-5 


75.1 


5.97-10 


-4 


262.1 


2.2 


GEM-6 


-17.4 


3.4- 10-5 


72.6 


9.7-10- 


-5 


233.2 


0.4 


GEM-7 


-14.9 


1.9- 10-* 


75.1 


1.1 - 10- 


-4 


196.4 


0.8 


GEM-8 


-14.9 


1.3- 10-^ 


75.1 


2.2- 10- 


-5 


251.6 


0.1 


GEM-9 


-14.9 


-5.7- 10-5 


75.1 


2.7-10- 


-4 


267.0 


1.0 


GEM-10 


-14.9 


2.6- 10-5 


75.1 


1.8-10- 


-4 


246.9 


0.7 


EGM96 


-14.9 


3.3-10-5 


75.1 


7.6-10- 


-5 


231.4 


0.3 



Table 7.b. The orientation of the ellipsoid of inertia 

As it is observed in Table 7.b, for the geopotential models considered, the 
longitude of Ox axis of the triaxial ellipsoid of inertia is about —15° and the 
longitude of Oy axis has a value close to 75°, except for the model GEM — 6 
for which ~ —17° and Xb — 73°. From the values obtained, except for the 
model GEM-5, it is determined the mean ellipsoid of inertia with the principal 
moments A = 8.010935639 - ■ , B = 8.011108377 - lO^^/cg - m^, C = 

8.037333747- lO^'^A:^ - and with the orientation = -15°. 2, As = 74°.8. 
The value of H obtained from the mean values of the principal moments of 
inertia, namely H — 0.00327369, is the same as in 7 . 

In Figure 1, the position of the inertial pole on the surface of the Earth is 
represented with respect to the conventional international pole Pq for the nine 
models of geopotential considered. Since the angle 7 is a small angle, then, in a 
Cartesian reference XPiyf in the tangent plan at Pqi with the axis P^X tangent 
to the Greenwich meridian, the pole of inertia Pi will have the polar coordinates 
(Ac, 7)- 

It is remarked that the pole of the ellipsoid SE-2 when 7 = coincides with 
the conventional international pole Pqj the other positions of the pole remaining 
in the neighborhood, except for the pole of inertia corresponding to GEM — 5. 
The coordinates of the mean pole P are Ac — 209°. 7 and 7 ~ 0".5. Therefore 
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the mean polar axis differs from the rotation axis by 0".5. The mean pole P 
deviates approximately by 15 meters from the pole of rotation. 





• 


. 1 " 




P 
O 




• 


1 1 i- 




2" 


1" 0",5 





Y 



'X, G 



Figure 1. The position of the inertial pole on the surface of the Earth 

The results obtained show that the approximation made in the paper |3] is 
satisfied and the improved values for the principal moments of inertia A, B, C 
are obtained. On the other hand, it is better emphasized the fact that the polar 
axis of inertia is located in the neighborhood of the Earth's rotation axis. For 
the geopotential models considered, the longitudes of the axes Ox and Oy of 
the triaxial ellipsoids of inertia have concordant values. 
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